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The statistics of triggered star formation: 

An overdensity of massive YSOs around Spitzer bubbles 
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ABSTRACT 

We present a detailed statistical study of massive star formation in the environment 
of 322 Spitzer mid-infrared bubbles by using the RMS survey for massive Young 
Stellar Objects (YSOs). Using a combination of simple surface density plots and a 
more sophisticated angular cross-correlation function analysis we show that there is a 
statistically significant overdensity of RMS YSOs towards the bubbles. There is a clear 
peak in the surface density and angular cross-correlation function of YSOs projected 
against the rim of the bubbles. By investigating the autocorrelation function of the 
RMS YSOs we show that this is not due to intrinsic clustering of the RMS YSO 
sample. RMS YSOs and Spitzer bubbles are essentially uncorrelated with each other 
beyond a normalised angular distance of two bubble radii. The bubbles associated with 
RMS YSOs tend to be both smaller and thinner than those that are not associated 
with YSOs. We interpret this tendency to be due to an age effect, with YSOs being 
preferentially found around smaller and younger bubbles. We find no evidence to 
suggest that the YSOs associated with the bubbles are any more luminous than the 
rest of the RMS YSO population, which suggests that the triggering process does not 
produce a top heavy luminosity function or initial mass function. We suggest that it 
is likely that the YSOs were triggered by the expansion of the bubbles and estimate 
that the fraction of massive stars in the Milky Way formed by this process could be 
between 14 and 30%. 
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1 INTRODUCTION 

Massive stars drive large amounts of energy into the sur- 
rounding interstellar medium (ISM) via their winds and 
ionising radiation, sculpting the ISM into a series of bub- 
ble and shell-like structures (Elmcgrccn 201 1|) . One partic- 
ular class of shells or bubbles are those surrounding ex- 
panding HII regions. The UV illumination excites emis- 
sion from Polycyclic Aromatic Hydrocarbons (PAHs) at the 
bubble rims and also heats dust grains within the HII re- 
gion. The former can be readily observed at 8 /im by IRAC 
onbo a rd th e Spitzer Space telescope (|Churchwell et all 
120061 . 120071 ) and the latter at longer wavelengths using 
e.g. M IPS on Spitzer jChurchwell etalj|2007l : I Watson et all 
120081) or PACS onbo ard the Herschel Space Observatory 
ijZavagno et al.ll2010al lbl). 

The expansion of these HII regions is of extreme in- 
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terest to studies of star formation as their expansion may 
trigger new generations of star formation into being within 
the molecular material surrounding the bubbles. There are 
two major triggering mechanisms that have been put for- 
ward so far: Radiative Driven Implosion and the Collect 
fc Collapse process. In Radiative Driven Implosion (RDI : 



i,ps 

Bertoldil 1 19891: iLefloch fe Lazarefj Il994l ; iMiao et al.l 120091 : 
Bisbas ct al.l |2009| , l201ll ) the expanding ionisation front 



of the HII region drives a D-type shock into molecular 
clouds surrounding the HII region, triggering the collapse 
of sub-critical clumps within the clouds. Theoretical mod- 
els of RDI can successfully explain the mo rphology of BRCs 
jMiao et al.ll2009l : lBisbas et al.ll2009l . l201ll ) and observations 
suggesting that star for mation is concentrated along the 
central axis of the c l ouds JSugitani. Tamura fe Qguralll999l : 
ISugitani et~aT1l2000l : iBisbas et al.ll201ll ). 

T he Collect fe Collapse process (Elmcgrccn fe Ladal 

1 19771 : IWhitworth et all 1 1994 IWhitworth fe Francis! |2002| ). 
on the other hand, does not require the presence of pre- 
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existing molecular structures. In this case the expansion of 
the HII region sweeps up the surrounding low density ma- 
terial into a shell surrounding the HII region. At a certain 
point (generally after a few Myr) this shell becomes self- 
gravitating, fragments and collapses to form dense molecular 
clumps that eventually collapse to form stars. The fragments 
formed in the Collect fc Collapse process tend to be mas- 
sive (a few hundred Mm . I Whitworth et al.|[l993) and so this 
process could naturally exp lain the hierarchica l nature of 
mass ive stellar clusters (e.g. iBastian et al.l [20051 ; lOev et all 
2005). It should be noted that neither mechanism excludes 
the other as a means of triggering star formation. Indeed, 
within a single HII region both RDI and Collect fc Col- 
lapse may operate together (|Deharveng. Zavagno fc Caplanl 
l2005h . 

However, recent models of the i onising feedback from 
massive stars (|Dale fe Bonnelll lioTlh suggest that neither 
mechanism is responsible and that the UV illumination from 
these stars simply erodes low density material rather than 
shaping its evolution and eventual collapse. In this scenario, 
the r ole of triggering in st ar formation is predicted to be min- 
imal. IWalch et al.l (|201lf ) find that the location of the dense 
clumps around the edge of HII regions (such as RCW 120) 
reflects the pre-existing cloud structure and their forma- 
tion does not require the Collect & Collaps e process. How- 
ever, r ather than neglecting triggering as in lDale fc Bonnelll 
l|201lf ). they suggest that stars may form by global implo- 
sion of the pre-existing structures (Enhancement of initial 
Density substructure and simultaneous Global Implosion, or 
EDGI). 

Observational studies of triggered star formation have 
so far mainly focused on photoionised globules (Bright- 
Rimmed Clouds or BR Cs, ISugitani. Fukui fc Qgural Il99ll ; 
ISugitani fc Qg ura 1994) found at the edges of optically visi- 
ble HII regions (iThompson et al.l2004| ; IUrouhart et al.l2004l; 



Morgan et all 120041 ; lUrquhart et all 120061 ; iMorgan et all 



2008) or on the rims of HII regions selected to show a rel- 



ative l y simple morphology dDeharveng. Zavagno fc Caplanl 
120051 ; IZavagno etail I2OO6I ; iPomares et al.l l2009h . Recent 
studies of the latter have been able to take advantage of 
the large catalogu e of infrared bubbles discov ered in the 
GLIMPSE survey l|Churchwell et al.ll200dl2007h . The prop- 
erties of these bubbles are consistent with expand ing HII re- 
gions jDeharveng et al.ll20ldl ; IWatson et al.ll2009l '). although 
perhaps more consistent with a ring ra ther than bubble mor- 
phology |Beaumont fc Williamsll2010l 'l. 

However, one difficulty with many of these studies is 
that they are often phenomenological in nature, concentrat- 
ing upon the visual identification of YSOs or protostars in 
regions where one might have a prior expectation that they 
may have been trigge red (e.g. IZavagno et al.l 120061 . 120071 ; 
iDeharveng et alj|2010h . Such studies cannot attack the cen- 
tral problem in triggered star formation, which is to iden- 
tify the origin of the discovered star formation (we refer to 
this as the "Origin Problem"). When trying to identify star 
formation as being triggered one must first exclude the pos- 
sibility that the star(s) would have formed spontaneously 
without the influence of the trigger. The Origin Problem is 
particularly intractable when considering individual objects 
— without a good understanding of the initial conditions 
involved it is almost impossible to categorise an individual 
star-forming region as being triggered or spontaneous. 



In this paper we attempt to move beyond the current 
phenomenological approach by carrying out a detailed sta- 
tistical stud^_ctf_^tar_jOTmation around the 322 Spitzer bub- 
bles of IChurchwell et all ^OOd ). While the Origin Problem 
is almost impossible to solve for individual star-forming re- 
gions, by considering the global properties of a large sample 
it may be possible to infer the presence of triggering in a 
statistical sense. These studies are made easier by the exis- 
tence of large, uniformly selected and well- understood star 
form ation surveys such as the RMS survey |Urouhart et al.l 
2008). The RMS survey has the goal of identifying every 
massive young stellar object (YSO) in the Milky Way and 
comprises an initial infrared selection followed by thorough 
multi-wavelength follow-up to rigorously classif y each ob- 
ject an d determine its physical properties (lUrquhart et al.l 



l2009al lbl; iMottram et~al]|2011al ; lUrquhart et al.ll201ll ). 



The RMS survey covers a greate r area than the 
GLIM PSE-I survey region in which the IChurchwell et ahl 
(2006) bubbles have been identified and is complete to YSOs 
of luminosity ^10 4 Lq out to the furthest bubble in the 
IChurchwell et ail (|2006T ) catalogue. Hence we can conduct a 
survey for recent massive star formation around all of the 
Spitzer bubbles (see Fig. [T] for an example). Note that the 
RMS survey does not include Galactic latitudes |Z| ^ 10° for 
reasons of co fusion and so we do not c onsider here the bubble 
catalogue of IChurchwell et al.l (|2007r i which is based solely 
upon the GLIMPSE-II survey region (i.e. \l\ 10°). On the 
other hand we must keep in mind the intrinsic subjective bi- 
ases in the catalogue of Spitzer bubbles, due to their manual 
by-eye identification by a number of independent observers 
ijChurchwell et al.ll2006l). These b i ases ar e discussed in detail 
in Section 2 of IChurchwell et al.l (|2006|). The current cata- 
logue is l ikely to be highly i ncom plete, particularly to small 
bubbles. IChurchwell et all (|2006h estimate a completeness 
on the order of ~50%, which is being borne out by early re- 
sults from the citizen science Milky Way Project (Simpson 
et al. 2012, in prep). 

Nevertheless, the IChurchwell et al.l l|2006T l catalogue of 
bubbles currently represents the most complete and well- 
studied catalogue of HII regions with a simple morphology 
that lends well to statistical studies of their YSO distribu- 
tion. In this paper we present such a study, paying partic- 
ular attention to the angular distribution of YSOs around 
the bubbles and potential differences within the YSO popu- 
lation. Our aims are to provide statistical evidence that the 
star formation associated with the bubbles may have been 
triggered and to investigate the properties of the global pop- 
ulation of bubbles and YSOs. In addition, the methods that 
we demonstrate in this paper will be readily applicable to 
the larger and more complete Milky Way Project bubble 
sample when it becomes available. 

The paper is structured as follows. In Sect. [5] we study 
the star forming environment of the Spitzer bubbles using a 
simple surface density approach, followed by a more sophis- 
ticated analysis of the angular cross-correlation function of 
RMS YSOs and Spitzer bubbles. Sect.|3]combines the results 
of the star forming environment analysis with the observed 
properties of the bubbles and YSOs to demonstrate that star 
formation is clearly enhanced toward the bubbles. We spec- 
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Figure 1. Thr ee colour GLIMPSE im age of a Spitzer bubble 
(N109 from the IChurchwell et alj|20o3 catalogue), showing the 
strong extended 8 fj,m PAH and 24 (im dust emission tracing the 
bubble rim. The colour coding in the image is 24 fim (red), 8 /im 
(green) and 4.5 p, m (blue). The 24 ^ m image is taken from the 
MIPSGAL sur vey iCarev et al.ll2q09h. N109 is one of the largest 
bubbles in the IChurchwell et alJ l|2006l ) catalogue, with a mean 
radius of 14.8', and is also the site of numerous smaller bubbles. 
The positions of two objects from the RMS YSO sample are in- 
dicated by green circles. 

ulate on the likely origin of this star formation and show 
that it is likely that the bubbles predate the formation of 
the RMS YSOs. Finally in Sect. [4] we present a summary of 
our conclusions and results. 



2 THE STAR FORMING ENVIRONMENT OF 
SPITZER BUBBLES 

2.1 The surface density of YSOs 

As a first approach to studying the distribution of RMS 
YSOs around Spitzer bubbles we simply measured the num- 
ber counts of RMS YSOs expressed as a function of angu- 
lar separation from the bubble centres. Our sample of RMS 
YSOs is comprised of the objects cla ssified as either YSO or 
UC HII in the RMS database^, see lUrquhart etafl <|2008h 
for a description of the RMS database and its classification 
system. YSO and UC HII sub-classifications both represent 
young, recently formed and predominantly massive stellar 
objects that enable us to trace the distribution of recent 
star formation around the bubbles. Hereafter we refer to 
this combined population as the RMS YSO sample. Within 
the a rea covered by the GLIMPSE I survey (|Beniamin et al.l 
120031 ) there are 846 ob jects within the RMS YS O sample and 
322 bubbles from the IChurchwell et al.l l)2006h catalogue. 

In order to account for the different angular radii of 
the bubbles we divided the angular separation of each RMS 
object from a particular bubble by the m ean radius of the 
bubb le ((R), column 9 in the catalogue of IChurchwell et alJ 
2006), i.e. expressing the angular separation in bubble radii 

2 http://www.ast.leeds.ac.uk/RMS 
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Angular offset (Normalised Bubble Radii) 

Figure 2. Histogram of the number counts of RMS YSOs (com- 
prising YSO and UC HII classifications from the RMS database) 
as a function of angular distance from the centre of Spitzer bub- 
bles. The distance is expressed in terms of normalised bubble 
radius. The number counts are scaled by the area of the annulus 
corresponding to each bin and thus represent a surface density of 
RMS YSOs. Error bars are determined via Poisson statistics. 

rather than arcminutes. Each bin represents an annulus 
around the centre of each bubble. The surface area of the 
annulus thus naturally increases with increasing radius and 
so to obtain the surface density of RMS YSOs we scale the 
counts in each bin by the surface area of each correspond- 
ing annulus. A histogram of these scaled number counts is 
shown in Fig. [2] 

Fig. [2] shows a clear peak in the number of RMS YSOs 
at a separation of 1 normalised bubble radius. At greater 
angular radii the number of RMS YSOs falls sharply, reach- 
ing a constant background level of ~ 3 sources by 2 bubble 
radii. Within an angular radius of 2 normalised bubble radii 
the number of RMS YSOs is demonstrably higher than at 
angular radii greater than 2 bubble radii. The surface den- 
sity of YSOs projected against Spitzer bubbles is thus higher 
than regions external to the bubbles, with a clear peak in 
the surface density projected against the the rims of the 
bubbles. 

The Spitzer bubbles are relatively elliptical, with typ- 
ical eccentricities between 0.6 and 0.7. As we normalise by 
the mean bubble radius (R) this will have the effect that we 
incorrectly calculate the true normalised radius of each RMS 
YSO from the bubble centre, potentially broadening the ob- 
served peak in surface density. The position angl es of the el- 
liptica l fits to the bubbles are not listed in IChurchwell et al.l 
(2006), but these measurements were kindly made available 
by Matt Povich (Povich, priv. comm.) so that we could ex- 
amine the effect of using the true radius of the bubble instead 
of the mean radius. We found that there is no significant dif- 
ference between scaling the distance of the RMS YSOs with 
the true bubble radius and the mean radius. This is more 
than likely due to the fact that the angular resolution of 
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our histogram in Fig. [2] is limited to 0.25 bubble radii by 
the need to obtain sufficient RMS YSOs in each bin. At this 
resolution the worst-case error in radius (i.e. between mean 
radius (R) and the semi-major axis a) is slightly larger than 
the width of one bin in Fig. [2] 

We further subdivided our RMS YSO sample into its 
constituent YSO and UC HII sub-samples to investigate 
trends in the separate distributions of YSOs and UC HIIs, 
for example in evolutionary status versus radius. We found 
that there is no significant difference between the two sub- 
samples, the histograms of separate YSO and UC HII sub- 
samples are indistinguishable from the combined RMS YSO 
sample. Again, this may be due to the limited sample statis- 
tics that we currently have, or this may indicate that gra- 
dients in evolutionary status around the bubbles are either 
not present or, if present, are on smaller angular scales than 
resolved by our study. 

In order to conf irm this result we also carried out the 
same analysis on the lRobitaille et al.l (|2008l ) catalogue of in- 
trinsically red sources selected from the Spitzer GLIMPSE 
survey. The Robitaille catalogue contains a much greater 
number of objects than the RMS survey, though at the ex- 
pense of_containhia^ii_by an uncertain fraction of AGB 
stars (Robitaille ct al. 2003). The surface density of objects 
from the Robitaille catalogue is shown in Figure [3] and dis- 
plays a similar distribution to the RMS YSO sample, with 
a higher surface density towards the bubbles that sharply 
drops off to a uniform background level. 

The background level is much higher than the RMS 
YSO samp le, as would be expected due to the higher surface 
density of iRobitaille et al.l (|2008l ) i ntrinsically red source s 
compared to the RMS catalogue ( Urauhart et all [2008). 
The distribution of IRobitaille et all ( 20081 ) objects does not 
peak at 1 bubble radius, but instead exhibits a relatively 
flat distribution out to 1 bubble radius. As the RMS cat- 
alogue is constrained to star forming objects at an early 
evolutionary sta t e (YS Os and UC HII regions), whereas the 
iRobitailie et al l i|200Sl ) catalo gue is not, this may indicate 
the presence of an evolution ary gradient acro s s the bubbles. 
Further classification of the IRobitaille et al.l (|2008r i sample 
and investigation of their star-forming nature are required 
to prove this hypothesis. 

We must also explore the possibility that there may 
be an intrinsic bias i n the distribution of both the RMS 
and IRobitaille et al.l (|2008l ) catalogues around the Spitzer 
bubbles due to the common mid-infrared bands used to 
detect both the bubbles and RMS/Robitaille objects. Al- 
though the bubbles are principally identified via their ex- 
tended PAH emission at 8 /im and the RMS YSOs and 
Robitaille intrinsically red sources are predominantly point 
infrared sources, the complex mid-infrared environments of 
the bubbles may lead to a bias in the identification of point 
sources at their rims. We investigate this possibility by ex- 
amining the distribution of 6.7 GHz methanol masers drawn 
from the Methanol MultiBeam (MMB) Survey (|Green et all 
2009) around the Spitzer bubbles. 6.7 GHz methanol masers 
are thought to exclusively tr ace young sites of massive star 
formation ('e.g. lMentenlll99ll '). and thus allow us to trace the 
distribution of massive YSOs around the bubbles indepen- 
dently of their mid-infrared emission. 

The MMB survey currently occupies a longitude range 
between / = 186 and / = 20, i.e. excluding the range 
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Figure 3. Histogram of the number counts of IRobitaille et alj 
(2008) Intrinsically Red Objects as a function of angular distance 
from the centre of Spitzer bubbles. The distance is expressed in 
terms of normalised bubble radius. The number counts are scaled 
by the area of the annulus corresponding to each bin and thus 
represent a surface density. Error bars are determined via Poisson 
statistics. 



20 ^ I ^ 186, and so only the bubbles in the south- 
ern GLIMPSE survey region are presently covered by the 
MMB survey. The individual masers in the MMB catalogue 
have had their positions interferometrically determined to 
sub-arcse cond precision and the maser det e ctions are re- 
>orted inlCaswelll (120091). [Green et al.l ll2009l).ICaswell et all 
2010l ). lGreen et al.l (|201oh . lCaswell et al.l l|201ll ) and Green 
et al. (2012, in press). We plot the surface density of 6.7 
GHz MMB masers around the southern Spitzer bubbles in 
Fig. HI 

Fig. [4] displays a v ery similar distribution of masers to 
that of RMS YSOs and IRobitaille et al.l i^POSf ) intrinsically 
red sources, albeit with larger error bars due to the smaller 
sample size. There is a clearly distinguished peak in the 
maser distribution at an angular offset of 1 bubble radius 
and the surface density of 6.7 GHz masers drops to a roughly 
constant background level beyond an offset of 2 bubble radii. 
The peak in the surface density of 6.7 GHz masers appears 
to be broader than the corresponding peak in the surface 
density of RMS YSOs, however the lower signal-to-noise of 
the MMB surface density makes it difficult to interpret this 
difference as a real effect. 

All three in d epend ently selected YSO catalogues (RMS, 
Robitaille et al.1 (|200Sf > red sources and MMB 6.7 GHz 
masers) display very similar surface density distributions 
and we thus conclude that the increase in surface density 
of YSOs towards the bubble rims is a real effect. Given 
this similar behaviour between catalogues, and the currently 
more compre hensive knowledge of t h e properties of the RMS 
YSOs (e.g. lUrquhart et al.l I2OIII ; iMottram etHl l2011al ; 
lUrauhart et al.1 l2009al lbf)~ we restrict our further analysis 
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Figure 4. Histogram of the number counts of MMB 6.7 GHz 
masers as a function of angular distance from the centre of Spitzer 
bubbles. The distance is expressed in terms of normalised bubble 
radius. The number counts are scaled by the area of the annulus 
corresponding to each bin and thus represent a surface density. 
Error bars are determined via Poisson statistics. 
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Figure 5. The angular cross-correlation of the RMS YSO sample 
and the catalogue of Spitzer bubbles as a function of normalised 
bubble radius. 



to the RMS YSO catalogue. Although the iRobitaille et all 
(2008) red source catalogue has greater sample statistics and 
is likely to be predominantly comprised of YSOs there is a 
much greater likelihood of contamination by AGB stars and 
other non-YSO types than in the RMS YSO catalogue. Sim- 
ilarly, the lower numbers of objects in the MMB 6.7 GHz 
maser catalogue favours the continuation of our study using 
the larger and much more studied RMS YSO sample. 

2.2 The angular cross-correlation of bubbles and 
YSOs 

As a refinement of our simple surface density ap- 
proach we also investigated the distribution of RMS 
YSOs around the bubbles using an angular two-point 
cross-correlation analysis, a technique more commonly 
used to determine the clustering properties of galaxies 
(e.g. ISmith Boyle fe Maddoxl Il995l; Ichirlanda et al.l l200d : 
IWang et al.l l201ll ; iBradshaw et alj 1201 lT ) . The correlation 
function defines the probability of finding a population of ob- 
jects at a particular angular separation from a different sec- 
ond population. Here, we use d the RMS YSO s a mple as our 
first population (-Di) and the lChurchwell eT al. (2006) cata- 
logue of infrared bubbles for our second population (D2). 
We calcul ated the angular cross- correlation using the es- 
timator of lLandv fc Szalavl |l993), modified for the cross- 
co rrelation betwee n popu lation 1 and 2 using the equation 
of IBradshaw et all (|201lh . i.e. 

IQ\ N D 1 D 2 — N Dl R 2 — N Rl D 2 + N Rl R 2 
W W - 7j (tj 

where Nd 1 d 2 represents the normalised number counts 
at an angular separation of 8 of RMS source-bubble pairs, 



Nd 1 r 2 and Nr 1 d 2 the counts of real and random catalogues 
of RMS source-bubble pairs (and vice- versa), and Nr 1 r 2 the 
counts between two random catalogues of RMS objects and 
Spitzer bubbles. As in Sect. 12.11 we scaled 8 to the radius 
of the individual Spitzer bubble in each pair. To avoid in- 
troducing high levels of noise through the randomly gener- 
ated catalogues we performed 50 realisations of each cat- 
alogue, taking the mean of the results to determine uj(8). 
The errors o n lj(8) were calculated by a bootstrapp i ng ap - 
proach (e.g. iGhirlanda et al.l 120061 ; IBradshaw" et al.ll201ll ). 
The Spitzer bubble catalogue was divided into 100 randomly 
chosen bootstrap catalogues each matching the original cata- 
logue size (with replacement). The angular cross-correlation 
was determined for each bootstrap catalogue and the result- 
ing la error in uj(8) is given by the standard deviation in 
u>(8) from the 100 random bootstrap samples. 

The resulting angular cross-correlation function is plot- 
ted in Figure [5] which reveals almost exactly the same dis- 
tribution as seen in the surface density distribution shown in 
Figures [2H1 The RMS YSO sample is found to be strongly 
correlated with the Spitzer bubble catalogue, particularly at 
the radius corresponding to the rim of the bubbles where the 
correlation peaks. This peak is significant at the 9er level. 
The cross-correlation drops sharply beyond the peak at 1 
bubble radius and beyond a distance of 2 bubble radii the 
cross correlation decreases to essentially zero. This indicates 
that the probability of finding an RMS YSO near a Spitzer 
bubble is markedly greater at an angular radius of 1 bubble 
radii, and that beyond an angular distance of 2 bubble radii 
the RMS YSO population are essentially uncorrelated with 
the presence of a Spitzer bubble. 



6 



2.3 The angular autocorrelation of RMS YSOs 

Finally we investigate the clustering within the RMS YSO 
sample, in order to determine whether the previous results 
in Sect. 12.11 and 12.21 are simply due to an intrinsic angu- 
lar clustering scale within the RMS catalogue that happens 
to correspond to the typical angular size of a Spitzer bub- 
ble. JKM_jnediar^_hubble_j^dius of all the bubbles listed in 
the IChurchwell et alj l|2006l ) catalogue is 1.1 arcmin, with 
a large observed range stretching between 0.14 arcmin and 
14 arcmin. However, few bubbles possess extreme radii and 
two thirds of the sample have radii between 1-3 arcmin. If 
the RMS YSO sample are naturally clustered on this typi- 
cal angular scale then the intrinsic clustering may mimic the 
apparent overdensity of YSOs observed in the surface den- 
sity histogram and angular cross-correlation plot (Figs. [2] & 

E0>. 

To compare the clustering of the control sample to that 
of the RMS YSO sample associated with bubbles we calcu- 
late the autocorrelation of the samp le as a function o f angu - 
lar distance, using the estimator of lLandv fc Szalavl (|l993l ) 



i-0 



N DD - 2N DR + N RR 

Nrr 



(2) 



Ndd, Nor and Nrr represent the normalised num- 
ber counts of data-data, data-random and random-random 
pairs respectively. We calculate the autocorrelation for three 
separate samples: the entire RMS YSO sample, those asso- 
ciated with Spitzer bubbles (i.e. lying within 2 bubble radii 
of a particular bubble), and those not associated with any 
Spitzer bubbles (i.e. lying more than three bubble radii from 
all bubbles). For simplicity we respectively refer to these 
samples as the full RMS YSO sample, the bubble-associated 
YSO sample and the control sample. 

We chose a value of 2 bubble radii for the radius of 
association due to the steep fall off in surface density and 
angular cross-correlation beyond this radius. A total of 116 
RMS YSOs are found within 2 bubble radii of a Spitzer 
bubble and 629 RMS YSOs are found at an angular dis- 
tance greater than 3 bubble radii from any bubble. As in 
Sect. 12.21 we constructed random catalogues of each of the 
three samples and performed 50 realisations of each random 
catalogue to avoid introducing higher levels of noise. 

The autocorrelations for these three samples of RMS 
YSOs are shown in Fig. El The behaviour of the full RMS 
YSO sample (solid dots in Fig. [6]) shows a classic peak to- 
wards smaller angular scales and a decrease towards larger 
angular scales. This implies that the full RMS YSO sample 
is strongly clustered on scales of ~1' or less. The bubble- 
associated YSO sample displays a markedly different be- 
haviour, being anticorrelated on all angular scales except for 
a small positive correlation at 2'. Finally, the control sample 
shows a correlation on small angular scales similar to the full 
sample, although much weaker and with a flatter fall-off to 
large angular scales than the full RMS YSO sample. Inter- 
estingly there is a minor peak in the autocorrelation function 
at 2', though this is not a statistically significant detection. 
At angular scales ^2' the autocorrelation functions of the 
full RMS YSO sample and the control sample are identical. 

Clearly the autocorrelation function of the bubble- 
associated RMS YSOs is very different to the other two sam- 
ples. This implies that the peaks in the YSO surface density 
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Figure 6. The angular autocorrelation function of RMS YSOs 
for the full sample of RMS YSOs (solid dots), a "control" 
sample of RMS YSOs that lie m ore than 3 bubble radii from 
all of the Churchwell et all l|200rJl bubbles (open squares), and 
the sample of RMS YSOs that lie within 2 bubble radii of 
the IChurchwell et alj l|200rj ) bubbles(open circles). Ict error bars 
are shown and are calculated using the bootstrap replacement 
method. 



and the YSO-bubble angular cross-correlation seen towards 
the bubble rims are not due to intrinsic clustering within 
the full RMS YSO sample. The YSOs in the full and control 
samples are much more highly correlated (i.e. clustered) on 
smaller angular scales than the bubble-associated YSO sam- 
ple. On the majority of angles the bubble-associated sample 
are anticorrelated, which implies that there is a shortfall in 
the numbers of YSOs associated with the bubbles at these 
angular separations compared to a random sample. The ex- 
ception to this is at an angular scale of 2', which corresponds 
closely with the median bubble diameter of 2.2'. This posi- 
tive correlation suggests that we may be seeing the signature 
of RMS YSOs located on either side of the bubble, and the 
anticorrelation implies that the YSOs are not found on an- 
gular scales smaller or larger than this. 

We do see a minor secondary peak at 2' in the au- 
tocorrelation functions of the full and control YSO sam- 
ples, which at first sight suggests that a fraction of these 
YSOs are correlated at the same angular scale as the me- 
dian bubble diameter. We expect this behaviour in the full 
YSO sample (which obviously includes YSOs associated 
wit h the bubbles). However , due to the incompleteness of 
the IChurchwell et al.l l|2006h catalogue caused by the man- 
ual search procedure that was used in its construction, we 
cannot exclude the presence of contaminating bubbles in our 
control sample. Thus the autocorrelations that we measure 
may be artificially enhanced for the control sample, for ex- 
ample at the secondary peak at 2'. However, this peak is 
not statistically significant and the full and control samples 
are clearly more strongly correlated at small scales than the 
bubble-associated sample. This suggests that the intrinsic 
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clustering of RMS YSOs does not produce the enhancement 
of YSOs projected against the rims of the Spitzer bubbles. 
Further work on the much more complete Milky Way Project 
bubble sample (Simpson et al. 2012, in prep) would aid this 
analysis. 



3 DISCUSSION 

In the following we combine our results on the YSO sur- 
face density distributions, the YSO-bubble angular cross- 
correlation and YSO angular autocorrelation together with 
the measured properties of the bubbles and RMS YSOs. Our 
particular aims are to investigate the star formation environ- 
ment of the bubbles in order to determine whether the YSO 
population is significantly enhanced near the bubbles and 
if there are any discernible differences in the population of 
YSOs found near the bubbles when compared to the entire 
sample. Ultimately we would like to place statistical con- 
straints on the star formation associated with the bubbles 
that can inform current and future models of triggered star 
formation. 

3.1 The properties of Spitzer bubbles associated 
with RMS YSOs. 

Taking the radius of association between YSO and bubbles 
to be 2 bubble radii we find that a total of 116 YSOs and UC 
HIIs fro m the RMS YSO sample are associated with bubbles 
from the lChurchwell et a l. (2006) ca talogu e. The number of 
bubbles from the lChurchwell et al.l (|2006h catalogue associ- 
ated with one or more RMS YSOs is 72, which corresponds 
to 22 ± 3% of the bubble catalogue (where the error is calcu- 
lated using Poisson statistics). Of the 72 bubbles associated 
with RMS YSOs, 30 bubbles are associated with more than 
one RMS YSO or UC HII, ranging from 2 to 5 objects per 
bubble (and with a mean of 2.6 objects per bubble for the 
multiple matches). 

The majority of the IChurchwell et all (|2006r i bubbles 
are not associated with RMS YSOs. This does not im- 
ply that the remaining bubbles are devoid of surrounding 
star formation, but merely that any star formation which 
is present is of sufficiently low luminosity to fall below the 
RMS detection thresh old (which is typical ly JSlOOO Lq at 
distances of a few kpc, Urquhart et al. 2011). Unfortunately 
the majority of the lChurchwell et al. I (|2006h bubbles do not 
have measured dist ances that are free from the kinematic 
distance ambiguity (Dc harveng et al"1l2010t ) and so we can- 
not determine the individual completeness limit for each 
bubble. Howeve r, the typical distances of the bubbles range 
from 2-13 kpc (|Deharveng et al.ll20ld ) and so we can con- 
fidently say that the present RMS study of bubbles is com- 
plete to YSOs with a lumino sity ^ 10 4 for th e most 
distant bubbles in the sample jUrquhart et alj|201ll ) and to 
YSOs with luminosity ^ 1000 L0 for bubbles located at the 
typical distance of a few kpc. 

As the luminosity of a B3 star is ~550 Lq 
l|Mevnet fc Mae dcr 2000) we thus identify massive star for- 
mation associated with the bubbles. Hence, 22 ± 3% of 
the bubbles are associated with massive star formation 
in the MYSO or UC HII re gion phase. This is a simi- 
lar fraction to that found by iDeharveng et all (|2010T l and 



IWatson. Hanspal fc Mengistul |2010l ) who respectively found 
18% and 20% of their bubble samples to be associated with 
either ultracompact HII regions, YSOs or 6.7 GHz methanol 
masers. The remaining 78 ± 3% of the bubbles may be as- 
sociated with low to intermediate mass star formation, but 
without a more sensitive survey we cannot confirm this hy- 
pothesis. The currentl y underway Herschel Hi-GAL survey 
of the Galactic Plane l|Molinari et al.|[2010al lbh will provide 
such a sensitive survey of the entire bubble sample and it 
would be advantageous to revisit the bubbles with the more 
sensitive Hi-GAL data when it is available. 

We searched for differences in the properties of the bub- 
bles that are associated with RMS YSOs and those that are 
not, in order to try and identify differences in the properties 
of bubbles that are associated with massive star formation 
and those that are not. The two most important measured 
properties of each bubble are the size (strictly the mean ra- 
dius of the bubbles) and thickness of the diffuse mid-infrared 
emission comprising the bubble. For both of these proper- 
ties we determined the mean values for the 22% of bubbles 
that are associated with at least one RMS YSO (i.e. within 
2 bubble radii) and the remaining 78% of bubbles that are 
not associated with an RMS YSO. 

The mean radius of bubbles that are associated with 
an RMS YSO is 3. 6' ±0.4', compared to the mean radius of 
unassociated bubbles of 4.6'±0.3'. The mean thickness of 
bubbles that are associated with RMS YSOs is 0.92' ±0.08', 
again compared to the unassociated bubble thickness which 
is 1.18' ±0.07'. Quoted errors are the standard error on the 
mean. In order to determine the significance of these dif- 
ferences in mean radius and thickness we performed a two 
sample unequal variance (heteroscedastic) t test on each pair 
of samples (RMS YSO-associated bubbles and unassociated 
bubbles). The t tests return probabilities of 1% and 0.8% 
respectively that the mean radius and thickness of YSO- 
associated bubbles and unassociated bubbles are drawn from 
populations with the same mean. 

Those bubbles that are associated with massive star 
formation thus tend to be both smaller and thinner than 
those bubbles that are not associated with massive star for- 
mation. These results are not significant at a level of 3cr or 
greater and, combined with the potential biases in the bub- 
ble population discussed earlier, should be interpreted with 
caution. However, it is instructive to speculate on what may 
be the physical cau ses behind these observed differences. 
I Weaver "et al. ( 1977) predict that the radius of a wind-blown 
bubble should increase much faster than the thickness of the 
swept up shell when the bubbles are in their first expansion 
stage. Hence, bubbles with small radii and thinner shells 
should be younger th an larger b u bbles with thicker shells, 
as also suggested by iDale et al.l {2009). Of course, we do 
not have physical distances for the majority of our sample 
and are thus dealing with angular radii rather than physical 
radii. Hence the bubbles that we have identified as small in 
angular size may just be the more distant members of the 
sample. However, over the whole sample of bubbles these 
effects should average out and our tentative results suggest 
that it is the younger bubbles within the sample that are 
more likely to be associated with massive star formation. 

Theoretical models of shell frag mentation (the 
collect and collapse process: IWhitworth et al" I Il994l ; 
IDale. Bonnell fc Whitworthl l2007f ) suggest that fragmenta- 



s 



tion of the swept up shell of a bubble tends to occur on 
timescales of one to a few Myr. Very little detailed study 
of the bubble lifetimes has currently been made. The few 
bubbles that have been stu died to date have dyna mical 
lifetimes of 0.5 to a few Myr (jWatson et al.ll200Sl , MMh and 
many of the sample are likely to be HII region s powered by 
late O to early B stars dDeharveng et al" 2010l; iBania et all 



Anderson et all l201lf lT 



l20ld ; iBeaumont fc Williams! I201C 
which again have main sequence lifetimes on the order of a 
few Myr to a few tens of Myr . So the picture that massive 
star formation tends to be associated with smaller and 
younger bubbles is largely consistent with the predictions 
of the collect & collapse models that star formation should 
happen on ~ Myr timescales. However, if these bubbles 
are shown to be much younger than 0.5 Myr then this 
would point towards the star formation being caused by 
the implosion of pre-existing density structures (e.g. the 
Enhancement of initial Density sub structure and simu lta- 
neous Global Implosion proposed bv lWalch et al .1 [20 1 lh . or 
perhaps sim ply pre-existing untrig gered star formation as 
suggested bv lDale fc Bonnelll (|201lh . 

In the following subsection we concentrate upon the sur- 
face density of YSOs around the bubbles, with the aim of 
showing that there is a statistically significant overdensity of 
YSOs and that it is unlikely that the majority of the YSOs 
formed spontaneously. 



3.2 An overdensity of YSOs around Spitzer 
bubbles 

It is clear from the surface density plots shown in Figures 
[2]H]that there is a significantly enhanced surface density of 
YSOs found at the rim of the bubbles. For the RMS YSO 
sample shown in Fig. [2] the YSO surface density is 19. 1± 
3.6 YSOs/unit area at an angular distance from the bubble 
centre of 1 bubble radius, some 4<r above the mean surface 
density of YSOs at angular distances greater than 2 bub- 
ble radii (3.1±0.2 YSOs/unit area). This result is confirmed 
by the angular cross-correlation function of the RMS YSO 
sample shown in Fig. [S] where the cross-correlation function 
peaks at an angular distance of 1 bubble radius with a sig- 
nificance of 9cr. In the immediate environment of a Spitzer 
bubble the highest probability location to find an RMS YSO 
is projected against the rim of the bubble. 

Moreover, it is clear from inspecting Fig. [2] that the 
surface density of YSOs is not only enhanced at an angular 
offset of 1 bubble radius, but that it is enhanced over the 
entire angular scale of the bubbles out to an angular offset of 
2 bubble radii. We can see this by comparing the mean sur- 
face density of YSOs "inside" 2 bubble radii and "outside" 
2 bubble radii. The mean surface density of YSOs within an 
angular offset of 2 bubble radii is 8.9±1.7 YSOs/unit area 
compared a value of 3.2±0.2 YSOs/unit area at an angular 
offset of 2 bubble radii or greater. A two sample unequal 
variance (heteroscedastic) t test of these two subsamples re- 
turns a probability of only 0.4% that these two subsamples 
are drawn from populations with the same mean. Hence we 
have demonstrated that there is a statistically significant 
overdensity of massive YSOs associated with the bubbles 
compared to the background, with an enhanced probabil- 
ity of finding these YSOs projected against the rim of the 
bubbles. 



What do these results imply? Firstly, there is a greater 
concentration of massive star formation towards the bub- 
bles than in the wider environment. This result is confirmed 
by the surface density of MMB 6.7 GHz masers (see Fig.Q, 
which trace a YSO population independently of mid- infrared 
emission. A greater concentration of star formation towards 
the bubbles implies that the bubbles are either efficient at 
producing YSOs, or that they are found in regions of high 
YSO surface density. This is the classic chicken and egg sce- 
nario applied to massive star formation: do the bubbles pre- 
cede the high surface density of YSOs, or does the high 
surface density of YSOs precede (or occur simultaneously 
with) the formation of the bubbles? 

Before considering this question more fully, we must 
bring in the second of our results - that there is an enhanced 
probability of finding YSOs projected against the rim of the 
bubbles (i.e. at an angular offset of 1 bubble radius). By in- 
specting the autocorrelation of the RMS YSOs we showed in 
Section [2]3] that this effect is not likely to be due to intrinsic 
clustering within the RMS sample on similar angular scales 
to the bubble radii. The ancillary question raised by this re- 
sult is: why are the YSOs more likely to be found projected 
against the rim of the bubbles, i.e. what is special about the 
bubble rims? 

The bubble rims are traced by 8fim PAH emission which 
originates from the photon-dominated region between the 
ionisation front being driven out by the HII region within 
the bubble and the surrounding neutral medium. The rim of 
the bubbles thus shows the interface between HII region and 
surrounding neutral gas. For a spherical bubble morphology 
one would expect the column density of gas to be greater 
at the bubble rims due to the greater path length through 
the neutral material towards the rims. So at first glance, the 
higher surface density of YSOs projected against the bubble 
rims may simply reflect the higher column density at the 
rim of the bubbles, i.e. the YSOs trace molecular column 
density. 

However, while the sample of bubbles that have been 
observed at relatively high angular r esolution in molecu- 
lar lines (|Beaumont fc Williams! |2010) do show a peaked 
molecular column density profile at a normalised bubble 
radius of 1, the column density falls off much less sharply 
than the YSO surface density . Inspection of Figure 2 from 
IBeaumont fc Williams! fgOlO) shows that at a normalised 
bubble radius of 1.5 the CO intensity can be roughly half 
of that at a normalised radius of 1. This suggests that the 
YSOs may not trace the column density distribution, al- 
though much closer scrutiny of the bubbles in a non-optically 
thick tracer is required to confirm this hypothesis. Moreover, 
the CO contrast betw een the centre of the bubbles and their 
rims is often extreme (| Beaumont fc Williams![2O10h whereas 
the YSO surface density within an angular offset of 2 bub- 
ble radii is everywhere higher than the background level. 
Thus we cannot confidently say that the YSO surface den- 
sity traces the gas column density around the bubbles. 

The YSO surface density is strongly peaked at an off- 
set of 1 bubble radius and decreases sharply beyond this 
value. Beyond an angular offset of 2 bubble radii the sur- 
face density of YSOs is essentially undistinguishable from 
the background level. The angular cross-correlation function 
shows a similar steep drop — beyond an angular distance of 
2 bubble radii the bubbles and RMS YSOs are essentially 
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uncorrelated. The implication of this is that whatever causes 
the rise in YSO surface density is closely related to the rim 
of the bubbles. The bubble radius is a dynamic value and 
expected to increase over time as stellar winds or radiation 
pressure causes the bubbles to expand. Combined with this 
is the fact that the massive YSOs and UC HII regions iden- 
tified by the RMS survey typically tend to have lifetimes 
around a few 10 4 to a few 10 5 years (jMottram et aljfeoilbl l 
and so should trace very recent star formation. 

The sum of these pieces of information leads us to con- 
clude that it is likely that the bubbles predate the YSOs. If 
the bubbles formed in an environment with a high surface 
density of YSOs (e.g. in the turbulent highly fragm ented 
initial conditions suggested bv lDale fc Bonnelj|201ll 'l. then 
the distribution should not peak at the rim of the bubble 
as the bubble rad ius is time-dependent. Similar arguments 
have been used bv lPreibisch et ajj (|201lf ) for YSOs detected 
at the edges of shells in Carina. Also in this case the extent 
of the enhanced YSO surface density should also not be re- 
lated to the current radius of the bubble - why are bubbles 
found in regions of enhanced YSO surface density occupy- 
ing twice their current angular radius? Finally, the relative 
timescales of the massive YSOs and those required for the 
expansion of the bubbles imply that the YSOs formed after 
the bubbles. We thus conclude that a significant fraction of 
the YSOs seen against the rim of the bubbles were likely 
triggered by the expansion of the bubble. 

A greater understanding of the dynamical timescales for 
the expansion of the bubbles and also the molecular envi- 
ronment of the bubbles are required to confirm this hypoth- 
esis. Pinpointing the YSO formation to have occurred after 
the bubble was formed is crucial to disentangling cause and 
effect in the star formation surrounding the bubbles. Cur- 
rently, only a few bubbles have had their dyna mical lifetimes 
estima ted and more studies similar to those of lWatson et alJ 
(2009) are required over a larger sample of bubbles. Com- 
paring the YSO distribution to the gas distribution is also 
crucial to investigate differences in the population of YSOs 
at the rims of Spitzer bubbles, for example to determine 
whether the star formation efficiency is enhanced at the 
bubble rims. Finally, it will also be instructive to apply 
the same statistical tools that we have used in this paper 
to th e latest generation of triggered star form ation models 
(e.g. IWalch et all l201ll : bale fc Bonnel l201ll ) in order to 
see whether the surface density of sink particles in the SPH 
simulations matches that of YSOs around observed bubbles. 



3.3 The luminosity function of RMS YSOs 
associated with bubbles 

The RMS survey has determined lumi nosities for their entire 
sample of YSOs and UC HII regions (|Mottram et al.ll2011al : 
lUrquhart et al.ll201ll ). and so in this subsection we seek to 
identify differences in the luminosity function between those 
YSOs that are associated with bubbles (i.e. within an angu- 
lar offset of 2 bubble radii) compared to the full population. 
As both luminosity functions are essentially power laws with 
turnovers caused by incompleteness at low luminosities, one 
must take care that any differences between the two sam- 
ples are not primarily due to differences in the completeness 
of each sample. To avoid this issue we cut both samples 
at a luminosity of I0 4 Lq, i.e. the minimum YSO luminos- 
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Figure 7. The luminosity function of the entire RMS YSO sam- 
ple (solid line) and those YSOs that lie within an angular offset 
of 2 bubble radii from a Spitzer bubble (dashed line). 
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Figure 8. The cumulative distribution function of the luminosi- 
ties of the entire RMS YSO sample (solid line) and those YSOs 
that lie within an angular offset of 2 bubble radii from a Spitzer 
bubble (dashed line). 



ity at which the RMS catalogue is complete ljUrauhart et al.l 
l201ll ). We plot these truncated luminosity functions in Fig. [7] 
and the corresponding cumulative distribution functions in 

Fig. m 

Inspecting Fig. [JJshows that the RMS YSOs associated 
with the bubbles do not tend to have a higher luminosity 
than the rest of the population. The two luminosity func- 
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tions are essentially identical. A Kolmogorov-Smirnov test 
on the two luminosity functions yields the result that there is 
a 27% probability that the two distributions are drawn from 
the same sample, and thus we cannot reject the null hypoth- 
esis with any significance. If the YSOs associated with the 
bubbles were triggered (as we argue in Sect. 13. 2p then we 
conclude that the triggering process does not result in stars 
with an appreciably different luminosity function. 

This result must be contrasted against the many 
inferences that have been made to date which suggest 
that the triggering process may produce stars with a 
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models for swept-up shells around massive stars predict a 
top-heavy mass distribution of fragment s, which in turn 
may l ead to a top-heav y IMF distribution (jWhitworth et all 
1 1994 iDale et al.l 12009 ). Observations of molecular clouds 
associated with HII regions also suggested that IRAS 
point sources associated with the clouds nearer the HII 
regions are more luminous tha n those found to be more 
distant <|Yamaguchi et al.l 1 1999 ). However these results are 
not consistent with the largely constant form of the IMF 
observed over Galactic scales. If the triggering process does 
induce a top-heavy I MF then the resu lts of triggering do not 
dominate the IMF l|Dale et al.ll2009l ). Here we have shown 
that if these YSOs are triggered, then the triggering process 
results in a luminosity function that is indistinguishable 
from the full YSO sample. 



3.4 The fraction of massive stars in the Milky 
Way that may have been triggered 

Estimating the impact of the triggering process in global 
Galactic star formation is a vital part of determining a sim- 
ple prescription for star formation that can be applied to 
galaxy evolution models. Because we have combined the re- 
sults of two large are a and relatively u n biase d surveys (the 
RMS survey and the IChurchwell et all (|2006l ) bubble cata- 
logue) we are in a position to try and estimate the contribu- 
tion of triggered star formation to the Galaxy's population 
of stars. Before d oing so, we must stress th e major caveat in- 
volved - that the lChurchwell et alj (|2006l ) bubble catalogue 
is likely to be highly incomplete. Thus, we do not attempt 
a detailed treatment at this stage and simply infer a lower 
limit to the fraction of massive stars in the Milky Way that 
may have been triggered. 

As discussed in Sect. 13.11 we find 116 YSOs and UC 
HIIs from the R MS YSO sample that are associated with 
bubbles from the IChurchwell et all (|2006l ) catalogue. The 
bubble-associated YSOs form 14% of the 846 objects con- 
tained in the RMS YSO sample over the GLIMPSE sur- 
vey region. The RMS survey is complete for essentially all 
massive YSOs with luminosities in excess of 10 4 Lq out to 
the furthest kinematic distance in the bubble sample (~14 
kpc). Thus, assuming that the incompleteness of the Church- 
well catalogue dominates over the fraction of the bubble- 
associated YSOs that were triggered, we estimate that at 
least 14% of the massive stars in t he Miilky Way cou l d hav e 
been triggered. If, as sugg ested by IChurchwell et al.l (|2006h . 



the Spitzer bubble catalogue is ~50% incomplete then the 
true fraction of triggered massive stars could be up to ~30%. 



4 SUMMARY AND CONCLUSIONS 

We have carried out a detailed statistical study of young 
stel lar objects (YSOs ) foun d nearby Spitzer bubbles from 
the IChurchwell et all (|2006l ) catalogue. Our main conclu- 
sions arc summarised below: 

(i) The su rface densities of YSO s and UC HIIs from the 
RMS survey ([Urauhart et al . 2011), intrinsically red sources 
from iRobitaille et al.l il2008l ) and MMB 6.7 GHz methanol 
masers ( Green et al.l 2009 | ) are enhanced towards Spitzer 
bubbles from the Churchwell et al. I (|2006l ) catalogue. The 
surface density of all three YSO catalogues peaks toward 
the projected angular radius of the bubbles, with a peak 
surface density of 4<r above the mean background level for 
RMS YSOs and UC HII regions. The mean surface density of 
YSOs associated with the bubbles is overdense with respect 
to the surrounding mean background at the 3<r level. 

(ii) The angular cross-correlation function of RMS YSOs 
and Spitzer bubbles shows a similar behaviour to the sur- 
face density, with a 9cr peak in the cross-correlation function 
at an angular offset of 1 bubble radius. In the immediate 
environment of a Spitzer bubble the highest probability lo- 
cation to find an RMS YSO is projected against the rim of 
the bubble. RMS YSOs and Spitzer bubbles are essentially 
uncorrelated beyond an angular distance of 2 bubble radii. 
Examination of the autocorrelation functions suggests that 
these effects are not caused by intrinsic clustering within the 
RMS YSO sample. 

(iii) Most Spitzer bubbles are not associated with 
massive YSOs: 22 ± 3% of the Churchwell cata- 
logue are associated with RMS YSOs. This fraction is 
consistent with smaller studies (|Deharveng et al.l l20ld ; 
I Watson. Hanspal fc Mengistullioiol ). 

(iv) Spitzer bubbles associated with RMS YSOs tend to 
possess both thinner rims and smaller angular radii than 
bubbles that are not associated with RMS YSOs. We in- 
terpret this tendency to be due to an age effect, with RMS 
YSOs forming around younger and smaller bubbles. 

(v) The different relative timescales for the formation of 
Spitzer bubbles and RMS YSOs, and the strong peak in sur- 
face density and cross-correlation at the rims of the bubbles 
lead us to conclude that a significant fraction of the RMS 
YSOs were triggered by the expansion of the bubble. 

(vi) We find no significant differences in the luminosity 
function of RMS YSOs associated with Spitzer bubbles com- 
pared to the entire RMS YSO population, which suggests 
that the triggering process does not result in a top heavy 
luminosity function or IMF. 

(vii) We estimate from the fraction of RMS YSOs associ- 
ated with bubbles and the incompleteness of the bubble cat- 
alogue that the lower limit for the fraction of massive stars 
in the Milky Way that could have been triggered is 14%. 
If the lChurchwell et al.l ((2006) catalo gue is 50% incomplete 
then the upper limit to the fraction of massive stars that 
could have been triggered may be up to ~30%. Therefore 
this mode of massive star formation ought not to be ignored 
when considering star formation on Galactic scales. 
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We must stress that these results are based on the 
current and largely incomplete sample of known Spitzer 
bubbles, which were identifi ed by manual searches of the 
GLIMPSE image database (IChurchwell et all 120061 ). How- 
ever the groundwork for a much more comprehensive cata- 
logue of bubbles is currently being laid by the Milky Way 
Project, which aims to identify many more bubbles by a 
systematic citizen science survey of the GLIMPSE images 
(Simpson et al 2012, in prep). In ad dition, the Herschel Hi- 
GAL survey l|Molinari et alj|2010bl) holds the promise of a 
much more complete YSO catalogue reaching fainter lumi- 
nosities than the RMS Survey. 

Statistical studies of the type that we have presented 
here will be of increasing importance in the age of large 
scale surveys of the Milky Way, and offer the prospect 
of being able to directly compare the observed distribu- 
tion of YSOs to the predicted distributi on of sink particles 
in triggered star form ation models (e.g. IWalch et alj|201ll ; 
iDale fc Bonnelll l201ll ) . With these advances the study of 
triggered star formation will finally move beyond the phe- 
nomenological stage to be able to make direct predictions. 
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